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ABSTRACT 

We derive improved system parameters for the HD 209458 system using a 
model that simultaneously fits both photometric transit and radial velocity ob- 
servations. The photometry consists of previous Hubble Space Telescope STIS 
and FGS observations, twelve /-band transits observed between 2001-2003 with 
the Mt. Laguna Observatory lm telescope, and six Stromgren b + y transits ob- 
served between 2001-2004 with two of the Automatic Photometric Telescopes 
at Fairborn Observatory. The radial velocities were derived from Keck/HIRES 
observations. The model properly treats the orbital dynamics of the system, 
and thus yields robust and physically self-consistent solutions. Our set of system 
parameters agrees with previously published results though with improved accu- 
racy. For example, applying robust limits on the stellar mass of 0.93-1.20 M , we 
find 1.26 < -Rpianet < 1-42 Rj up and 0.59 < M p i anct < 0.70 M Jup . We can reduce 
the uncertainty on these estimates by including a stellar mass-radius relation 
constraint, yielding _R p i anct =1.35±0.07 Rj up and M planet = 0.66 ±0.04 M Jup . Our 
results verify that the planetary radius is 10-20% larger than predicted by irradi- 
ated planet evolution models, confirming the need for an additional mechanism to 
slow the evolutionary contraction of the planet. A revised ephemeris is derived, 
T = 2452854.82545 + 3.52474554E" (HJD), which now contains an uncertainty in 
the period of 0.016 s and should facilitate future searches for planetary satellites 
and other bodies in the HD 209458 system. 

Subject headings: stars: individual, HD 209458 — extrasolar planets 

1. Introduction 

Of more than 150 recently discovered extrasolar planets 1 , HD 209458b is the first known 
to transit its star (Henry et al. 2000; Charbonneau et al. 2000). Combined with the spec- 
troscopic radial- velocity curve, photometric observations of transits allow high-precision de- 
termination of system parameters such as the inclination, planetary radius and mass. For 



1 A current tally can be found at the California and Carnegie Planet Search home- 
page, http://exoplanets.org, or at the J. Schneider Extrasolar Planet Encyclopaedia, 
http: / /www. obspm.fr / encycl/ encycl.html. 
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example, using Hubble Space Telescope photometry of HD 209458b, Brown et al. (2001) 
derived estimates of the inclination i = 86°.6 ± 0°.14 and planetary radius -R p i an et =1-347 
±0.060 Rj up ; with the additional estimate that the host star's mass is M star =1.06 ±0.13 M , 
Cody & Sasselov (2002) find M planct =0.69 ±0.02 M Jup . 

With a measured radius and mass, HD 209458b observations are extremely important 
in refining the theory of irradiated extrasolar giant planets (EGPs). Most notably, the 
observed radius of HD 209458b is ~10-20% larger than models predict (Guillot & Showman 
2002; Bodenheimer, Lin, & Mardling 2001; Bodenheimer, Laughlin, & Lin 2003; Baraffe 
et al. 2003; Laughlin et al. 2005). Guillot & Showman (2002) use evolutionary models of 
irradiated EGPs to argue for an additional heat source acting on HD 209458b. They found 
that the radius of HD 209458b can be produced with the transport of ~1% of incident stellar 
flux into the lower atmosphere as kinetic energy (i.e. winds). As HD 209458b is probably 
tidally locked, one side is perpetually illuminated, and high-speed winds would be expected 
to transfer heat from the day side to the night side (Showman & Guillot 2002; Cho et al. 
2003; Menou et al. 2003). Bodenheimer, Lin, & Mardling (2001) suggest that continuing 
tidal circularization of HD 209458b is heating the planet, inflating it to a radius larger than 
predicted for a circular orbit. However, a non-circular orbit is strongly disfavored by recent 
timings of the secondary eclipse of HD 209458b by Deming et al. (2005). Observations of the 
transiting exoplanet TrES-1 (Alonso et al. 2004; Sozzetti et al. 2004; Laughlin et al. 2005), 
indicating that its radius is consistent with models that do not invoke additional heat sources, 
support the suggestion of Burrows et al. (2004) that the large radius of HD 209458b may 
be anomalous. The radii estimates of the OGLE-discovered exoplanets are also consistent 
with models that include only irradiation. Determining the planetary radius of HD 209458b 
more accurately would assist in resolving this discrepancy. 

An extremely precise ephemeris can allow the inference of additional bodies in the 
system by their gravitational effects on the transit times and slight asymmetries in the 
transit light curve (Sartoretti & Schneider 1999; Agol et al. 2005; Holman & Murray 2005). 
The transit of HD 209458b affords the opportunity to make highly accurate and precise 
measurements of the planet's orbital period. The ability to predict and observe HD 209458b 
transits at 1-second precision would place strong constraints on the masses of satellites or 
additional planets in the system. Brown et al. (2001) point out that an Earth-mass satellite 
of HD 209458b would alter the time of mid-transit by up to 13 seconds. An ephemeris precise 
at the 1-second level would facilitate tests for such timing displacements. The precision of 
the STIS transit timings obtained by Brown et al. (2001) allowed the exclusion of satellites 
larger than 3 Earth masses at the 3a level by analysis of the timing displacements of the four 
transits they observed. Even if no moons are present around HD 209458b, the technique 
of transit timing can be readily applied to future transiting systems. Some extrasolar giant 
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planets reside in the habitable zones of their parent stars, a region where the star's insolation 
is such that liquid water could exist on the surface of a planet. Terrestrial-size moons of 
giant planets in such an orbit could in theory be habitable (Williams, Kasting, & Wade 
1997). While the mid-transit timing perturbation caused by a hypothetical satellite orbiting 
HD 209458b is on the order of 10 s, the perturbation grows as the cube root of the planet's 
orbital period, so for a planet in the habitable zone the perturbation due to a satellite could 
be tens of seconds to minutes in duration. The only currently available methods of detecting 
such satellites require transits; refining modelling techniques using HD 209458b thus builds 
a solid foundation for characterization of future discoveries of transiting extrasolar planets, 
e.g. the Kepler mission. The most widely used ephemeris by Robichon & Arenou (2000), 
based on Hipparcos photometric data and the radial- velocity curve (Mazeh et al. 2000), 
contained uncertainties in the period of ±1.21 seconds. Propagating to the present epoch 
yields an uncertainty in the time of mid-transit of over ±16 minutes. There exists, therefore, 
a need to revise the ephemeris, as the predicted time of mid-transit continuously accumulates 
errors in period. Even with a period determination accurate to 1 second, the large number 
of cycles (~100 per year) results in nearly 2 minutes of uncertainty accumulating after only 
a year. 

In this work, we fit a single, self-consistent model to observations of twenty-seven 
HD 209458b transit events in four bandpasses, and to more than three years of high-precision 
radial velocities. In §2 we briefly describe the five data sets used in this study, §3 outlines 
the modelling procedure, and in §4 we present and discuss the system parameters. 

2. Observations 

2.1. Mount Laguna Observatory (MLO) 

Twelve /-band transits of HD 209458b were observed with the 1-meter telescope at 
Mount Laguna Observatory (MLO). A log of these and all other photometric data is pre- 
sented in Table 1. Observing runs in 2002 and 2003 consisted of a 4 to 6 hour series of 
1-second exposures on a Loral 2048 x 2048 CCD. A field of view of 5x6.5 arcmin was used, 
which included three faint comparison stars. Observations in 2001 used 4-second exposures 
and a field of view of 5.9 x 5.9 arcmin that only included TYC 1688-1903-1 as a comparison 
star, and so this star was used as the comparison on all nights. 

The brightness of HD 209458 (V=7.645) required defocusing the telescope slightly to 
avoid saturating the CCD; as such, large photometry apertures (~3x FWHM or 9 arcsec) 
were employed. Light curves were obtained using standard differential photometry with one 
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comparison star. The photometry was then binned by a factor of 5-10 to a median time 
resolution of 177 s to help reduce the scintillation noise of the short exposures, yielding a 
total of 1149 points. The error bars on the binned points were calculated in the following 
way: uncertainties based on propagation of the IRAF-generated error estimates were used 
to compute the inverse-variance weighted mean per bin. Then this error bar was boosted by 
a/ (xl) if the reduced xl °f the scatter of the points about the mean in the bin was greater 
than I. 

For most transits, slight tilts (~0.5%) in the light curves were evident, likely due to 
systematic calibration errors. Comparison star variability was ruled out, as the tilts did 
not occur in all nights. Color-dependent differential extinction was also excluded, as the 
correlation between the tilts and the airmass was not consistent. These tilts were corrected 
by masking the transit and fitting a line to the out-of-transit light curve. This linear fit was 
then subtracted from the data, rectifying the light curves. Finally, under the assumption 
that the out-of-transit light curve is constant, the root mean square (RMS) deviation of 
the out-of-transit observations was compared with the mean uncertainty of these data. We 
found the error bars to be slightly underestimated, almost certainly due to systematic errors 
unaccounted for in the data reduction. We boosted all the uncertainties by 11% to make the 
uncertainties consistent with the out-of-transit RMS deviation. 



2.2. Hubble Space Telescope STIS 

Brown et al. (2001) obtained and analyzed extremely high-precision observations of 
four HD 209458b transits using the Hubble Space Telescope Space Telescope Imaging Spec- 
trograph (STIS). A total of 20 Hubble Space Telescope orbits in four visits spanning 18 days 
yielded a total of 684 spectra, covering a wavelength range of ~582-638 nm. Binned over 
wavelength, the spectra yield photometry with a relative precision of about 1.1 x 10~ 4 per 
60 s integration. Following Brown et al. (2001) we omitted the first orbit of each of the four 
HST visits in our modelling because of a systematic 0.25% deficit in flux, leaving 556 points 
in the light curve. As with the MLO observations, the RMS deviation of the out-of-transit 
light curve about a constant was used to check for the accuracy of the error estimates. We 
boosted the uncertainties by 4.5% to make the error estimates agree with the RMS deviation 
out of transit. 
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2.3. Hubble Space Telescope FGS 

Five transits were observed by Schultz et al. (2003) using photomultiplier tubes (PMTs) 
in the Hubble Space Telescope Fine Guidance Sensor (FGS). The FGS, normally used for 
astrometry and pointing control, was for the first time used as a high-speed photometer on 
a bright star. The four PMTs in FGSlr were used with a sampling rate of 40 Hz, and a 
S/N~80 per 0.025 s exposure was achieved. The F550W filter was used, giving a central 
wavelength of approximately 550 nm. A time dependency in the FGS response was corrected 
by fitting a 5th order Chebyshev polynomial to the out-of-transit data for each transit and 
each PMT. The data were then placed into 80 s bins for consistency with the STIS data of 
Brown et al. (2001). Schultz et al. (2004) have used these FGS data and the STIS data of 
Brown et al. (2001) to obtain a precise set of system parameters: i?pi anet =1.367±0.043 Rj up , 
-R s tar=l-154±0.036 R , and inclination i=86°.525±0°.054, assuming a stellar mass of M star = 
1.1±0.1 M from Mazeh et al. (2000). An orbital ephemeris of T =2452223.895819±0. 000031 
HJD and P=3.52474408±0. 00000029 days was also obtained; the uncertainty in the period 
is only 0.025 s. The precision of these data allowed the exclusion of transiting satellites of 
HD 209458b down to 2 Earth radii. Such an object would cause a ~0.1% dip in the transit 
light curve, which was not seen in the FGS data. We trimmed a total of 19 data points 
from the FGS observations, as these points were > 4<r outliers from an initial model fit to 
the FGS data set. These points occurred at the beginning and end of each orbit, where the 
reliability of the FGS data is the poorest, probably due to HST "breathing." A total of 268 
FGS points were used in the modelling. The uncertainty estimates of the original Schultz et 
al. (2004) data were boosted by 38% to account for systematic errors and be consistent with 
the out-of-transit RMS deviation about a constant flux. 



2.4. Automatic Photometric Telescope (APT) 

Six transits were obtained in 2001-2004 with the T8 and T10 Automatic Photometric 
Telescopes (APTs) located at Fairborn Observatory in the Patagonia Mountains of southern 
Arizona (Henry 1999; Eaton, Henry, & Fekel 2003). The transits of HD 209458b were co- 
discovered in 2000 using data from the T8 APT (Henry et al. 2000). Precision photometers 
use dichroic mirrors to split the incoming light into two beams and two EMI 9124QB bi-alkali 
photomultiplier tubes to measure Stromgren b and y simultaneously. Differential magnitudes 
from the two passbands were combined into a single (b + y)/2 band for greater precision, 
which is typically 0.0012 mag for a single measurement. The comparison star was HD 210074 
(V=5.74 mag, F2V). These APT data had a time resolution of 95 s and provide a total of 
1426 observations. The uncertainties were boosted by 7% so as to be consistent with the 
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out-of-transit RMS deviation about a constant. 

2.5. Keck/HIRES Radial Velocities 

A set of 51 published (Henry et al. 2000) and unpublished Keck/HIRES radial velocity 
measurements was included in the model fits and is listed as Table 5 in the Appendix. These 
data were obtained with the HIRES echelle spectrometer at a resolution of R ~ 80000. An 
iodine absorption cell was used for wavelength calibration (Marcy & Butler 1992; Valenti, 
Butler, & Marcy 1995; Butler et al. 1996). The current data set includes 42 measurements not 
included in Henry et al. (2000), and covers a time span of 3.4 years. The mean uncertainty 
of these velocity measurements is 4.7 m s -1 . 

3. Modelling of Transit Photometry and Spectroscopy 

3.1. Outline of the Physical Model 

The Eclipsing Light Curve (ELC) code (Orosz & Hauschildt 2000) was used to model 
HD 209458 by simultaneously fitting all transit light curves and the radial velocities. ELC 
explicitly includes the orbital and rotational dynamics of the system, yielding robust and 
physically self-consistent solutions. In ELC, the surfaces of the star and planet are defined by 
equipotential surfaces in the rotating binary frame. By using Roche equipotential surfaces 
to determine the radii, ELC can account for any non-sphericity of the bodies. In practice, 
for HD 209458 the difference between the polar and equatorial radii is only ~0.007%, and 
can be neglected. ELC includes the effect of gravity darkening, though again because of the 
near sphericity of the star, this effect results in a negligible (~0.002%) change in temperature 
between stellar pole and equator. 

In the original version of ELC, two parameters called the "filling factors" (/star and 
/planet) were used to define the equipotential surfaces. The filling factors are defined as the 
ratios of radii to Roche lobe radii, such that / < 1 for detached systems and / = 1 for 
a Roche lobe-filling body. For given masses, the filling factors thus define the stellar and 
planetary radii. Unfortunately, this way of specifying the equipotential surfaces is far from 
optimal in the case of HD 209458. Since the two bodies are nearly perfect spheres, a change 
in the mass ratio Q = M star /M planet requires a change in / star and / p i an et to get the same 
transit profile. Thus it proved to be much more computationally efficient to have as the two 
parameters the fractional radius of the star R stSur /a and the ratio of the radii -R s tar/-R P ianet 
since the transit light curve is independent of the mass ratio when the equipotential surfaces 
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are defined this way. 

Once the surfaces of the star and planet are defined they are then divided into a grid of 
surface elements ("tiles"). For HD 209458 a very fine grid of surface elements was required: 
even with 400 latitude x 600 longitude tiles on the star we found numerical noise in the light 
curves of the order of 1 part in 10 5 , which is not good enough to model the precision of the 
Hubble Space Telescope light curves. To mitigate this effect, we employed a Monte Carlo 
subsampling of the partially eclipsed tiles, using 1000 random subsamples per tile. This 
technique allowed a much smaller number of grid elements and we adopted a final 250 x 240 
grid on the star and 80 x 80 grid on the planet. 

The intensity at each tile is approximated as a blackbody, accounting for limb darkening 
and correcting the effective temperature T e g for gravity darkening (von Zeipel 1924; Claret 
2000). Limb darkening was treated using the 2-parameter logarithmic prescription of van 
Hamme (1993): 

1(h) = 7o(l — x + xji — yiAn/i) (1) 

where h is the foreshortening angle of the grid element and x and y are the two limb-darkening 
coefficients for each bandpass. Blackbody intensities were employed for convenience, but also 
because the Kurucz model atmosphere tables are presently too coarse at the T e fr=6000 K 
temperature appropriate for HD 209458 (Mazeh et al. (2000), Cody & Sasselov (2002)). After 
correcting the intensity of each tile for limb and gravity darkening, the binary is then "turned 
in space" by a user-specified phase step. A phase step size of 0°.05 was chosen, corresponding 
to about 42.3 s, to match the high time resolution (80 s) and precision (~ 1.1 x 10~ 4 ) of the 
Hubble Space Telescope data. If a tile on the star is completely eclipsed by the planet, that 
hidden tile is not included in the calculation of total flux from the star. Partially eclipsed 
tiles are accounted for via the Monte Carlo method mentioned above. In this way, a model 
light curve is generated for the input parameters. 

In the original version of ELC, the scale of the binary was fully specified by using the 
inclination i, the orbital separation a, and the mass ratio Q. In the case of HD 209458, we 
found it more convenient to use the mass of the star M star and the projected semiamplitude 
of the star's radial velocity curve K star . For a given inclination i, the separation a and the 
mass ratio can be determined from M star and -K" s tar- Then given the orbital period P, the 
radial velocity of the star at each orbital phase is computed by summing the velocities of each 
stellar tile not eclipsed by the planet, again weighting for limb darkening and foreshortening 
(Wilson & Sofia 1976). In this way we fit the radial velocity observations, and do not just 
constrain the stellar motion to equal some adopted projected orbital velocity K star . 

Since ELC computes the projected radial velocities at each orbital phase, the Rossiter 
effect can be fit with our model. The Rossiter effect is a radial-velocity distortion which 
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occurs as a result of the planet blocking the approaching and receding limbs of the star during 
transit (Rossiter 1924). The Rossiter effect was first observed in HD 209458 by Bundy & 
Marcy (2000) and Queloz et al. (2000), and was used to deduce that the planet orbits in 
the same direction as the star's rotation. The amplitude of the Rossiter effect is directly 
proportional to the star's rotational velocity V rot sini and the size of the eclipsing object. 
As portions of the rotating stellar surface are blocked by the planet, the absorption line 
profiles are skewed in shape. Thus, the velocities reported here during the transit represent 
an apparent Doppler shift caused by the skewing of the original absorption lines, rather 
than a change in velocity. Because of this, we caution that the radial velocities through the 
transit as computed by ELC may not be identical to the radial velocities as measured in the 
stellar spectra via the cross-correlation method. A thorough investigation would require the 
convolution of a high signal-to-noise out-of-transit spectrum (or model spectrum) with a set 
of phase-dependent line-broadening functions generated during the transit. These synthetic 
spectra should then be cross-correlated with a template spectrum and the resulting simulated 
radial velocities compared with the ELC radial velocities. Such an investigation is beyond 
the scope of this paper, but we note that the Rossiter effect observations are in fact very 
well matched by the ELC model as seen in §4.2. 

A circular (e = 0) planetary orbit was assumed, consistent with a tidal circularization 
time of order 10 8 yr for HD 209458b (Bodenheimer, Laughlin, & Lin 2003) and radial velocity 
observations (Mazeh et al. 2000). Recent Spitzer Space Telescope observations by Deming et 
al. (2005) show the secondary eclipse occurs at orbital phase 0.5000 ± 0.0015, demanding a 
very low eccentricity. The planet is also assumed to be tidally locked (rotation synchronous 
with its orbit) so the rotational frequency parameter f2 p i anet , defined as the ratio of rotational 
frequency to the orbital frequency, was set to unity. The effective temperature of the planet 
was set to 1300 K, based on theoretical estimates of the night-side temperature of irradiated 
extrasolar giant planets (Guillot & Showman (2002); Deming et al. (2005) measure a bright- 
ness temperature of 1130 ± 150 K). At this temperature the flux ratio of the planet to star 
is approximately 3 x 10~ n in the V band, and so the planet is essentially invisible in our 
optical photometry. 

Finally, note that ELC computes the actual orbital path of the planet about the barycen- 
ter rather than simply a chord across the stellar disk. This includes the reflex motion of the 
star during transit, which for a 3.07 hr transit, amounts to a distance equal to 0.9% the ra- 
dius of the planet. Though these effects are subtle, if improperly treated they can be sources 
of systematic error and can bias the system parameters. By including the dynamics of the 
orbits, ELC should in principle yield more accurate, as well as precise, system parameters. 
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3.2. 



Observational Constraints 



HD 209458 is a single-lined binary system, and it is well-known that the radial velocity 
curve of the star sets the minimum mass of the planet: 



Once the mass function /(M p i anet ) is known, one needs the inclination and the mass of the 
star (or the ratio of masses) in order to find the mass of the planet M p i anet . Once the 
component masses are known, the orbital separation a can be found from Kepler's third law. 
Since the transit light curves constrain the ratio of the radii (-R s tar/-Rpianet) and the fractional 
radii (Rstm/a and -Rpianct/a), the radii of the star and planet in physical units are found once 
a is known. 

From a strictly observational point of view, the inclination, relative radii, and limb 
darkening are constrained quite well from the multi-color transit light curves. However, 
neither the stellar mass or the mass ratio is constrained by the transit or radial velocity 
observations and so the scale of the system is indeterminate (see also the discussion in Deeg, 
Garrido & Claret (2001)). Since we are interested in the physical size and mass of the 
planet, an additional constraint is required to break the degeneracy. Fortunately we do have 
constraints on the radius and mass of the host star, by matching stellar evolution models 
with the star's luminosity and temperature. 

Cody & Sasselov (2002) derive a stellar radius of 1.18±0.10 R using the apparent V 
magnitude, the Hipparcos parallax (Perryman et al. 1997), and a well-determined bolometric 
correction. The ~ 8% error on this stellar radius estimate yields an 8% error on the orbital 
separation and hence will yield an uncertainty of ±8% in the planetary radius. We include 
this stellar radius estimate as an observation to be matched by the ELC model. The model 
is directed toward this value by use of a straightforward \ 2 penalty for deviations from the 
estimate. In a similar fashion, we include the observation of the projected stellar rotation 
velocity V TOt sin i of Queloz et al. (2000) as a datum to be matched. Thus in addition to the 
transit photometry and radial velocities, we include two more observables, with the radius 
of HD 209458 "steered" toward 1.18±0.10 R©, and toward a V rot smi of 3.75±1.25 km s _1 . 
However, we emphasize that these are not constraints — the model is permitted to adopt 
values outside of these ranges, but at a cost in X '■ 

Cody & Sasselov (2002) also determined a robust limit on the stellar mass of 1.06±0.13 
M ; this uncertainty range includes the observational errors in temperature, luminosity, and 
metallicity as well as systematic errors in convection mixing-length and helium abundance. 
The ±12% uncertainty on the stellar mass translates to roughly a ±4% error on the planet 



/(M planet ) = 




2ttG " (M planct + M star ) 2 ' 
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radius. However, this range of allowable stellar mass is not quite a 1-a uncertainty, but more 
of a limit on the range of possible mass. For this reason, and also because the stellar mass 
determination is more model-dependent than the other parameters, we do not include the 
stellar mass as an observed parameter. 

The relative uncertainty in the orbital velocity K stai also maps into an uncertainty in 
orbital separation and hence radii. The correspondence is 1:1, so an uncertainty of 2.1% (see 
§4.2) in K star gives a minimum uncertainty of 2.1% in radii. This 2.1% is the uncertainty in 
the radii if the mass of the star were exactly known. 



3.3. Fitting Procedure 

Using the above assumptions and conditions, all four photometric time series were mod- 
elled simultaneously with the radial velocities in order to obtain system parameters for 
HD 209458. The criteria defining the goodness of fit was the usual x 2 ■ The ELC model 
contains 15 free parameters: the orbital period of the planet P, the time of mid-transit To, 
the fractional radius of the star i? sta r/«, the ratio of the radii -R s tar/-Rpianet, the inclination, 
the semiamplitude of the star's radial velocity curve -ft'star, the ratio of rotational frequency 
of the star to orbital frequency, and eight wavelength-dependent limb-darkening coefficients 
(two per bandpass). Because of the degeneracies in the solutions, the stellar mass was not 
a free parameter in the models. Rather, we fix the stellar mass and optimize all other pa- 
rameters, then repeat with a different stellar mass. In this way we ran 69 models, stepping 
through stellar masses between 0.57 and 1.72 M . The list of other fixed input parameters 
and their values is given in Table 2. 

The weighting of the five data sets was determined by their uncertainties; no adjustments 
to the relative weights were made. In order to assure complete exploration of parameter 
space, a genetic algorithm based on the pikaia routine given in Charbonneau (1995) was 
used to find the global x 2 minimum (see Orosz et al. (2002) for details), and a downhill 
simplex "amoeba" (Press et al. 1996) was used to further examine the minimum. A simple 
grid search was then used to step through the stellar mass. 

For the purposes of this investigation, the input light curves were scaled in the following 
way: APT data were treated as Johnson B band light curves, FGS data were considered as 

V band, STIS data were considered R band, and MLO data were treated as / band. Using 
the Hipparcos V magnitude of 7.645 (Perryman et al. 1997) and the observed V — I, R — I, 
and B — V color indices (H0g et al. 2000), the out-of-transit fluxes were scaled to B = 8.18, 

V = 7.645, R = 7.287, and / = 6.985 mag. Although the effective bandpasses of the four light 
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curves are not equivalent to Johnson filters, these approximations are acceptable since ELC 
compensates by adjusting the limb darkening parameters. Also note that the bandpasses 
are treated independently, i.e. the relative fluxes between them are not constrained by the 
blackbody temperatures, so the out-of-transit scaling can in fact be arbitrary. 

4. Results and Discussion 

4.1. Light Curve Fits 

Figure 1 shows plots of all photometric data phase-folded and overlaid with the ELC 
model light curves. For this figure, the APT and MLO light curves were binned by a factor 
of 5 for clarity. The effect of color- dependent limb darkening is evident in the shape of the 
transits: the I-band light curve is considerably flatter during mid-transit than the B-band. 
The reduced x 2 of each data set are the following: MLO: 1.96, STIS 2 : 0.95, FGS: 1.27, and 
APT: 1.01. The reduced x 2 °f the entire fit to all the observations is 1.36 with 3434 degrees 
of freedom. The reduced x 2 values for the STIS and APT observations are excellent; for the 
MLO and FGS observations the x 2 is high, probably as a consequence of systematic errors 
in the data rather than deficiencies in the model. Normalized residuals ((O — C)/a) of the 
model fits are shown in Fig. 2 and, with the exception of the MLO data, show no pattern 
of correlated residuals. The MLO residuals do exhibit structure: just prior to mid-transit 
the data are generally above the model, and after mid-transit the data are below the model. 
We are confident this is an observational calibration problem, perhaps a residual remaining 
after the linear rectification. 



4.2. Radial- Velocity Fit 

Shown in Figure 3, the model radial velocities match the observations throughout the 
orbit. No obvious patterns indicative of a nonzero eccentricity are evident in the residuals. 
The reduced x 2 of this model fit to the radial velocities is 2.16, though removal of one outlier 
at phase 0.05 lowers the reduced x 2 to 1.80. The RMS scatter of the observations about the 
fit is 6.6 m s -1 , not much larger than the mean uncertainty of 4.7 m s _1 in the data. We 
derive a stellar reflex velocity -K" sta r = 82.7 ± 1.3 m s _1 (note that this is independent of the 



2 If we fit to the full STIS light curve, including the out-of-transit data obtained during the first HST 
orbit of the visit, the STIS reduced \ 2 is 1-54; this is a consequence of the systematically lower flux for these 
unreliable data, as described in Brown et al. (2001). 
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assumed stellar mass). This is in good agreement with the results of Henry et al. (2000) 
81.5±5.5 m s -1 , and marginally smaller than the values of Mazeh et al. (2000) 85.9±2.0 m 
s- 1 and Naef (2004) 85.1±1.0 m s -1 . 

Figure 3 also shows a detailed plot of the fit to the Rossiter effect, showing the excellent 
agreement between the ELC model and the observations. During the transit the RMS of 
the residuals is 3.9 m s _1 , whereas the out-of-transit RMS is 8.5 m s _1 . The data taken 
during transit were all obtained on the same night (JD 2451755), whereas the remaining 
data come from widely separated epochs spanning more than three years. Systematics such 
as the long-term stability of the spectrograph are thus most likely the cause of the increased 
scatter outside of transit. 



4.3. The Orbital Ephemeris of HD 209458b 

While extremely precise, the STIS observations only span 6 nearly consecutive HD 209458b 
orbital cycles and therefore are not ideal for determining the ephemeris. The FGS observa- 
tions span 1.3 years and hence place much tighter constraints on the period. By combining 
the STIS and FGS high-precision data with the ground-based observations spanning many 
cycles, very tight limits can be placed on the orbital period. We obtain a revised ephemeris 
of T = 2452854.82545 + 3.52474554E HJD (Table 3). The uncertainty in period is now only 
0.016 s, compared to 1.21 s (Robichon & Arenou 2000), an improvement by over a factor of 
70. A comparison of our ephemeris with previously published results is shown in the O-C 
diagram of Fig. 4, which displays the residuals of times of mid-transit from our ephemeris. 
We include most previously published mid-transit times (Charbonneau et al. 2000; Jha et 
al. 2000; Mazeh et al. 2000; Deeg, Garrido & Claret 2001; Schultz et al. 2003, 2004). Since 
we have included the STIS and FGS data in our fits, the estimates of T by Schultz et al. 
(2004, 2003) and Brown et al. (2001) are not independent and so we cannot use the O — C 
values in Fig. 4 to further improve the ephemeris. 

Our results are in general agreement with previous work, but some discrepancies exist; 
for example, our period differs from Schultz et al. (2004) by ^8 a (~ 0.126 sec). The majority 
of the uncertainty in the ephemeris lies in T ; the uncertainty in T is 750 times larger than 
that in the period. Clearly, if T can be determined as accurately as P, attaining 1-second 
precision in transit timings, and hence searches for satellites and other gravitationally per- 
turbing bodies, would be far less challenging. Despite their very high photometric precision, 
neither the FGS or STIS can obtain an uninterrupted observation of the 185 min transit due 
to HST's low Earth orbit. This makes determining the center of transit much more difficult, 
and hence more uncertain. What is needed, then, are observations of complete transits at 
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a precision comparable to the Hubble Space Telescope data, or a large set of lower-precision 
mid-transit times spanning many cycles. The Canadian MOST satellite, expected to achieve 
micromagnitude photometric precision, is well-suited for this task, as its orbit enables it to 
continuously observe HD 209458 for many transits (Rucinski et al. 2003) with exceptionally 
high precision. 

4.4. HD 209458 System Parameters 

Given only the transit light curves and radial velocities, the absolute physical scale of 
the system cannot be determined and the solutions are degenerate. In Fig. 5 we show the 
observed relationships between the planetary and stellar radii versus the stellar mass. The 
solid dark curves are the ELC solutions that best match the transit and radial velocity data, 

1 /3 

and as expected, the curve follows R oc M st .' ar . Formally, a fit to the model points yields the 
relation # planet = 1.310 x (M star ) a334 R Jup and R star = 1.112 x (M star ) - 333 R , and these 
are shown as the dark curves in Fig. 5. The planetary and stellar radii are observationally 
constrained to lie along these curves, though their positions on the curve are only weakly 
constrained. By including Cody & Sasselov's (2002) observationally-derived stellar radius 
of 1.18±0.10 R Q we break the degeneracy, but the uncertainty in the system scale remains 
large, giving a minimum uncertainty of ±8% in the planetary radius and 25% in the stellar 
mass. Our best fit value for the stellar mass using only the stellar radius constraint (along 
with the transit photometry, radial velocities and V ro tsini) is 1.25±0.35 M . 

If we apply the Cody & Sasselov (2002) stellar mass limit constraint, 1.06±0.13 M , 
we limit the solutions to the light gray region in Fig. 5. The smaller, darker gray region 
corresponds to the mass range given by Santos et al. (2004), 1.15±0.05 M Q . The vertical 
width of the gray regions correspond to the ±2.1% uncertainty in radius derived from the 
uncertainty in K star . With these mass range estimates, we can robustly bracket the planetary 
radius to the range 1.28 - 1.42 Rj up . Note that this is not a statistical 1-er confidence interval; 
solutions outside this range are highly disfavored. For comparison, we show the planetary 
radius estimates of Cody & Sasselov (2002), Deeg, Garrido & Claret (2001), Brown et al. 

(2001) , and Schultz et al. (2004). The latter three studies assumed a stellar mass of 1.10 
M Q to get the planetary radius, and so their radius ranges should be compared to the height 
of the gray region on our radius vs. mass curve. The box and diamond symbols show the 
locations of the points corresponding to the best stellar mass estimates of Cody & Sasselov 

(2002) and Santos et al. (2004). 

The lower panel of Fig. 5 shows the family of stellar radius vs. stellar mass solutions for 
HD 209458. This relationship comes solely from the observations of the transit and radial 
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velocities, and does not include any stellar astrophysics. From the transit and radial velocity 
observations alone, there is no difference in goodness of fit anywhere on this curve. The mass- 
radius relationship for HD 209458 given by Cody & Sasselov (2002), based on stellar models, 
is shown as the dot-dash line (bracketed by the uncertainty in the stellar radius). As Cody 
& Sasselov (2002) noted, this M-R relationship is nearly orthogonal to the transit curve, 
and thus introduces a strong constraint on the stellar mass and radius. The intersection of 
the Cody & Sasselov (2002) M-R relationship with our transit-derived relationship occurs 
at 1.09 M Q and defines our best estimate of the stellar mass. Incorporating the uncertainty 
in the mass-radius relation then yields these values: M star =1.093 ±0.092 M , i? star =1.145 
±0.056 R , and i? p i anct =1.350 ±0.066 Rj up , respectively, though again we caution that the 
uncertainties are not statistical 1-a confidence intervals. 

We estimate a slightly higher stellar mass and smaller stellar radius than Cody & Sas- 
selov (2002). This requires a slightly different set of adopted system parameters, e.g., we 
find the inclination to be 86.67 versus 86.1 degrees. But more importantly, our smaller stel- 
lar radius requires a higher temperature to produce the observed luminosity. Interestingly, 
Santos et al. (2004) estimate the temperature to be 6117±26 K, significantly hotter than the 
Cody & Sasselov (2002) or Mazeh et al. (2000) value of 6000±50 K. Ribas et al. (2003) also 
find a somewhat higher temperature, 6088±56 K, as do Fischer & Valenti (2005), 6099±44 
K. Even a small temperature increase of 88 K requires a 3% change in radius to maintain 
the same bolometric luminosity, and would reduce the radius from the Cody & Sasselov 
(2002) value of 1.18 R to 1.145 R ; a 117 K increase in effective temperature would drop 
the radius to 1.136 R ; both of these agree well with our stellar radius estimate, suggesting 
stellar temperatures greater than 6000 K are favored. These higher temperatures imply a 
slightly earlier spectral type, F9, than the usual adopted value of GO (see e.g. Gray, Napier 
& Winkler (2001) and Gray, Graham & Hoyt (2001)). 

A higher temperature would also mean a younger age, as can be seen in the stellar 
evolution tracks shown in Fig. 1 of Cody & Sasselov (2002): the observation box would 
be shifted to the left, and closer to the main sequence. However, the slightly younger age 
makes little difference to the planetary radius predicted from planetary evolution, since the 
shrinking of the planet vs. time is quite slow after 10 9 years e.g. see Burrows et al. (2001). 

As a check, we estimated the stellar radius from the Hipparcos parallax-derived distance, 
the observed magnitudes (assuming zero reddening), the observed surface gravity (Santos et 
al. 2004), an assumed effective temperature, and synthetic photometry computed from the 
NextGen stellar models (Hauschildt et al. 1999). These stellar radius estimates are shown 
in Fig. 5, for T cff =6000 K and 6117 K. Also shown are the stellar radius estimates of Cody & 
Sasselov (2002), Allende Prieto & Lambert (1999), and Ribas et al. (2003); all are consistent 
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within their uncertainties and with our "best" mass and radius estimate. In particular, the 
optical-IR synthetic photometry technique employed by Ribas et al. (2003) gives a stellar 
radius of 1.145±0.049 R in excellent agreement with our favored estimate. 

It is interesting that the Cody & Sasselov (2002) planetary radius estimate is notably 
larger than our value, despite our heavy reliance on their stellar radius estimate to set the 
physical scale of the system. Using our transit-derived radius vs. mass relation, their stellar 
radius estimate (1.18 R ) is inconsistent with their preferred stellar mass (1.06 M Q ); a 
larger stellar mass, 1.20 M , is required for consistency with our transit modelling. For 
completeness, we note that the radius vs. mass relationship derived by Deeg, Garrido & 
Claret (2001), R star = 0.34M star + 0.825 (±0.06), gives a stellar radius larger by ~0.05 R 
than our results, and hence their system scale is larger by approximately this amount. Their 
estimate was made prior to the exquisite Hubble Space Telescope transit observations and the 
tighter constraints those observations provide; our stellar radius vs. mass relation is similar 
to their lower 1-a limit on the relation. 

Just as with the planetary radius, the planetary mass depends on the adopted stellar 
mass. The transit and radial velocity data constrain the planetary mass to lie along a curve; 
formally we find M planct =0.620 x (M star ) a670 M Jup from our model fits. This is exactly what is 
expected for circular Keplerian orbits and mass ratio Q >> 1: M p i mict =(P/2nG) 1 ^ 3 (K star / sin?) x 
(M sta r) 2//3 - Applying the stellar mass-radius constraint as above, we find M p i anct =0.658 
±0.036 M Jup . 

In Table 4, we present the system parameters as a function of stellar mass for three 
representative cases, M star = 0.93, 1.09, and 1.19 M . These cases bracket the full range 
in stellar mass that is deemed acceptable, following Cody & Sasselov (2002). The overall 
reduced x 2 °f the fit to the 3434 observations (transit photometry, radial velocities, -R s tan 
and Vrotsini) was 1.364. Because of the degeneracy of the solutions discussed above, the x 2 
does not vary significantly across this mass range. 

Our system parameter values are generally in good agreement with those derived in 
previous work. Our estimate of the planetary radius is strongly constrained to be between 
1.26 - 1.42 Rj up , using the stellar mass limits of Cody & Sasselov (2002) and Santos et 
al. (2004). Applying the stellar mass-radius relation of Cody & Sasselov (2002), we get 
the limits: 1.35 ±0.07 Rj up . This radius is consistent with previously published results 
(Brown et al. 2001; Cody & Sasselov 2002; Schultz et al. 2004; Deeg, Garrido & Claret 
2001; Laughlin et al. 2005), and hence the discrepancy between the observed and theoretical 
models of extrasolar giant planets remains. The planetary radius is 10-20% larger than 
evolutionary models which only include irradiation (Chabrier et al. 2004); these give a radius 
of ~1.1 Rju P , equal to the bottom edge of the top panel in Fig. 5, and firmly ruled out. 
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However, Burrows et al. (2003) point out that the observed transit radius (radius where 
the planet's optical depth to starlight transmitted through the atmosphere along our line 
of sight is unity) is not the same as the theoretical 1 bar radius, and the difference can 
account for as much as 0.1 Rj up in the radius. When this effect is included, irradiated models 
can be marginally consistent with the extreme lower limit on the observed radius. As an 
alternative explanation for the planet's large radius, Bodenheimer, Lin, & Mardling (2001) 
suggested that eccentricity pumping by an undetected additional planet could provide the 
energy needed to slow the evolutionary contraction of HD 209458b. However, the radial- 
velocity residuals currently do not support a significant nonzero eccentricity. The recent 
observation of the secondary eclipse of HD 209458b by Deming et al. (2005) also strongly 
supports an eccentricity indistinguishable from zero, as the timing of the secondary eclipse is 
not significantly displaced from phase 0.5. Thus the observed radius of HD 209458b remains 
unsatisfactorily explained and continues to demonstrate deficiencies in our understanding 
of irradiated extrasolar giant planet evolution. Finally, with a mass of 0.66 ±0.04 Mj up and 
the above radius, the mean density of HD 209458b is 0.33±0.05 g cm~ 3 , about half the 
density of Saturn and one-fourth the density of Jupiter; HD 209458b is a very low density 
gas giant planet. 



5. Summary 

We have used the ELC code (Orosz & Hauschildt 2000) to determine the system param- 
eters of HD 209458 by simultaneously fitting the transit light curves and radial velocities. 
The observations consist of the Hubble Space Telescope STIS (Brown et al. 2001) and FGS 
(Schultz et al. 2004) light curves, plus 18 transits obtained over 4 years with the facilities 
at MLO and APT, along with Keck HIRES spectroscopic radial velocities. Our new esti- 
mates of the system parameters are generally in agreement with previous results, e.g., 1.26 
< -Rpianet <1.42 Rj up and 0.59 < M p i anct <0.70 Mj up . We stress that this range includes the 
uncertainty in the stellar mass. By applying the mass-radius relation of Cody & Sasselov 
(2002), we reduce the uncertainty: i? p i anct =1.35 ±0.07 Rj up , M planet =0.66 ±0.04 M Jup . Our 
results confirm that the planetary radius remains significantly larger (10-20%) than pre- 
dicted by irradiated planet evolution models e.g. see Chabrier et al. (2004). For the stellar 
parameters, we find M star =1.09 ±0.09 M and R star =l.l5 ±0.06 R . 

We have also obtained an orbital ephemeris with a period determination good to 0.016 s, 
over 70 times more precise than the period by Robichon & Arenou (2000). An ephemeris 
of this precision should facilitate future searches for additional bodies in the HD 209458 
system. By using ELC, a fully self-consistent dynamical model that includes subtle physical 
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effects not contained in previous work, we have reduced systematic errors. As observational 
precision increases, and detection methods become sensitive to smaller planets at larger 
orbital distances (ultimately to Earth analogs), the need for such exactitude in the modelling 
is warranted. 
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Table 1. Observation Log 



Instrument 


UT Date 


Start Time (HJD-2450000) 


M.LU 


onm i livi on on 
ZUUl Jun Zy-oU 


onon o a n 
zUoy.o4U 


M.LU 


zUUl Oct lo 


nine; a oa 
zlyo.ooo 




onno Tim i q 
ZUUz Jun lo 


z4oo.oob 


MLO 


onno A . . c 

ZUUz Aug 5 


o a m to 

2491.(00 


MLO 


onno A 1 o 

ZUUz Aug 12 


o a no 70/ 1 

z49o.7zb 


MLU 


onno o „-»-. 07 
ZUUz bep z7 


Z044.DZD 


MLO 


onno /~\ _ a. ,i 

ZUUz Oct 4 


2551.64U 


MLO 


onno /~\ „ a. 1 i 

ZUUz Oct 11 


o et cr o i 

zooo.olo 


i\,tt r\ 
MLU 


onno {~\ 1 c 
ZUUz Oct lo 


ZOOO.Ool 


MLO 


onno t~\ in 

ZUUz Dec 1U 


2618.585 


MLU 


onnu T.ii 07 
ZUUo Jul z< 


OO/IV r 7'7 a 

284/.770 


MLO 


2003 Aug 3 


2854.718 


STIS 


2000 Apr 25 


1659.744 


STIS 


2000 Apr 28-29 


1663.297 


STIS 


2000 May 5-6 


1670.336 


STIS 


2000 May 12-13 


1677.376 


FGS 


2001 Jun 11 


2072.274 


FGS 


2001 Sep 11-12 


2163.902 


FGS 


2001 Nov 10 


2223.793 


FGS 


2002 Jan 16 


2290.783 


FGS 


2002 Sep 30 


2548.085 


T8 APT 


2001 Oct 6 


2188.592 


T8 APT 


2001 Oct 13 


2195.579 


T10 APT 


2002 Oct 4 


2551.589 


T10 APT 


2002 Oct 11 


2558.581 


T10 APT 


2004 Sep 15 


3263.606 


T10 APT 


2004 Sep 22 


3270.597 



Table 2. ELC Model Input Constraints 



Parameter 


Value 


Reference or Reason 


T e g, star 


6000 K 


Cody & Sasselov (2002) 


T e ff, planet 


1300 K 


Guillot & Showman (2002) 


Orbital eccentricity e 


0.00 


Mazeh et al. (2000) 


Planet rotation/orbital frequency (fi p i ane t) 


1.00 


Assume tidal locking 


Radius of star i? s tar 

Rotational velocity of star V ro t sin i 


1.18±0.10 R 
3.75±1.25 km s" 1 


Cody & Sasselov (2002) 
Queloz et al. (2000) 
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Table 3. HD 209458b Ephemeris 



Parameter 


Value & Uncertainty 


T (HJD) 


2452854.82545 ± 1.35 x 1CT 4 


Period (days) 


3.52474554 ± 1.8 x 10" 7 



Table 4. HD 209458 System Parameters 



Parameter 


Estimate 






M st ar (Mq) 


0.93 


1.09 


1.19 


-Rstar (Rq) 


1.085 


1.144 


1.178 


M p i anet (M Jup ) 


0.593 


0.657 


0.697 


-^planet (R-Jup) 


1.279 


1.349 


1.388 


Inclination (degrees) 


86.668 


86.668 


86.668 


Kstar (m/s) 


83.0 


82.7 


82.7 


Orbital separation (AU) 


0.044 


0.047 


0.048 


Limb-darkening coefficients 








x(B) 


0.877 


0.874 


0.872 


y{B) 


0.285 


0.278 


0.276 


x(V) 


0.724 


0.724 


0.724 


y(v) 


0.327 


0.328 


0.328 


x(R) 


0.774 


0.775 


0.775 


y(R) 


0.434 


0.435 


0.436 


x(I) 


0.819 


0.818 


0.817 


y(i) 


0.637 


0.634 


0.632 


X 2 (3434 degrees of freedom) 


4682.3 


4681.4 


4681.4 


xl 


1.364 


1.363 


1.363 




Orbital Phase 



1 

Orbital Phase 



Fig. 1. — Light curves and model fit, in normalized flux units (with arbitrary offset) folded on 
our ephemeris. The orbital phases shown here span 6.7 hours and each tick mark corresponds 
to 25 minutes. Left panel: HST FGS and STIS observations. Right panel: Ground-based 
APT and MLO observations. 
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Fig. 2. — Left panel: Normalized residuals (x = (O — C)/a) of ELC model fit to FGS and 
STIS data sets. Right panel: Normalized residuals for APT and MLO data sets. 
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Fig. 3. — Top panel: Radial velocity curve and model fit. Bottom panel: Expanded view of 
the Rossiter effect during transit. [Editors: Please rotate figures clockwise 90 degrees.] 



-27- 



-i 1 1 r 



~ i 1 1 1 1 1 1 1 r 



o - 



1 5 



j_8 
■ 7 



i- 



m 

u -200 - 
i 

o 



-400 



__U 

-400 



I , , , , I , , , , I 

-300 -200 -100 

Cycle Number 



Fig. 4. — O-C diagram of published transit times. 1: Charbonneau et al. (2000), 2: Mazeh 
et al. (2000), 3: Jha et al. (2000), 4: STIS (Brown et al. 2001), 5: STIS (Schultz et al. 2003), 
6: estimated from Deeg, Garrido & Claret (2001) Figure 1, 7: FGS & STIS (Schultz et al. 
2004), 8: FGS (Schultz et al. 2003). Our result is shown at cycle 0. 
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Fig. 5. — Planetary and stellar radii versus stellar mass, as determined by the transit and radial 
velocity data. The shaded regions show acceptable solutions when a stellar mass constraint is 
used (light gray = Cody & Sasselov (2002), dark gray = Santos et al. (2004)). For comparison, 
several published radius estimates are shown. The dot-dash line in the lower panel is the mass- 
radius relationship of Cody & Sasselov (2002). Being nearly orthogonal to the transit-derived 
relationship, this strongly constrains the solutions. 
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Table 5. Appendix A. HD 209458 Radial Velocities 



JD-2440000 Velocity (m/s) Uncertainty (m/s) 



11341.120 


27.4 


4.3 


11368.941 


-14.9 


4.7 


11372.134 


-68.9 


5.4 


11373.056 


63.4 


4.2 


11374.055 


57.3 


3.8 


11410.012 


-46.9 


4.1 


11410.963 


-58.2 


4.3 


11411.933 


67.8 


4.7 


11438.808 


-75.1 


4.5 


11543.689 


-4.4 


5.8 


11543.693 


7.2 


6.2 


11550.691 


-13.1 


5.8 


11550.695 


-21.9 


5.7 


11551.696 


-69.8 


6.4 


11551.701 


-61.8 


5.9 


11552.703 


59.6 


6.3 


11679.107 


-30.1 


4.4 


11703.121 


-72.9 


4.1 


11704.098 


17.7 


3.8 


11705.105 


81.4 


4.4 


11706.102 


-41.1 


4.7 


11707.108 


-49.5 


4.1 


11754.975 


39.0 


4.3 


11755.019 


29.6 


3.7 


11755.025 


22.4 


4.6 


11755.041 


25.1 


3.8 


11755.053 


56.8 


4.2 


11755.068 


55.2 


3.9 


11755.074 


48.7 


3.8 


11755.079 


42.5 


3.5 


11755.085 


37.3 


4.1 


11755.090 


33.4 


3.8 


11755.096 


18.2 


3.8 


11755.101 


10.4 


4.2 


11755.107 


0.0 


3.9 


11755.113 


-8.4 


4.2 


11755.118 


-19.3 


4.2 


11755.124 


-28.9 


4.2 


11755.129 


-30.6 


4.0 


11755.972 


-71.5 


4.5 


11792.791 


73.2 


5.0 


11882.707 


-71.4 


6.3 


11883.720 


21.5 


6.6 


11899.731 


-4.7 


5.4 


11900.721 


-62.5 


5.7 
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Table 5 — Continued 



JD-2440000 


Velocity (m/s) 


Uncertainty (m/s) 


12063.112 


-32.4 


4.7 


12102.004 


-37.7 


4.9 


12446.128 


-37.0 


4.5 


12514.965 


42.7 


4.1 


12535.845 


14.8 


5.0 


12575.730 


97.8 


6.7 



